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ABSTRACT 

We make detailed theoretical predictions for the assembly properties of the Local 
Group (LG) in the standard ACDM cosmological model. We use three cosmological 
N-body dark matter simulations from the Constrained Local Universe Simulations 
(CLUES) project, which are designed to reproduce the main dynamical features of 
the matter distribution down to the scale of a few Mpc around the LG. Additionally, 
we use the results of an unconstrained simulation with a sixty times larger volume 
to calibrate the influence of cosmic variance. We characterize the Mass Aggregation 
History (MAH) for each halo by three characteristic times, the formation, assembly 
and last major merger times. A major merger is defined by a minimal mass ratio of 
10: 1. 

We find that the three LGs share a similar MAH with formation and last major 
merger epochs placed on average « 10—12 Gyr ago. Between 12% and 17% of the 
halos in the mass range 5 x lO^^/i~^M0 < Mh < 5 x lO^^/i~^M0 have a similar MAH. 
In a set of pairs of halos within the same mass range, a fraction of 1% to 3% share 
similar formation properties as both halos in the simulated LG. An unsolved question 
posed by our results is the dynamical origin of the MAH of the LGs. The isolation 
criteria commonly used to define LG-like halos in unconstrained simulations do not 
narrow down the halo population into a set with quiet MAHs, nor does a further 
constraint to reside in a low density environment. 

The quiet MAH of the LGs provides a favorable environment for the formation of 
disk galaxies like the Milky Way and M31. The timing for the beginning of the last 
major merger in the Milky Way dark matter halo matches with the gas rich merger 
origin for the thick component in the galactic disk. Our results support the view that 
the specific large and mid scale environment around the Local Group play a critical 
role in shaping its MAH and hence its baryonic structure at present. 

Key words: galaxies: haloes; cosmology: theory; methods: N-body simulations 



1 INTRODUCTION 

Observations of the Milky Way (MW) and the galaxy M31 
shape to a great extent our understanding of galaxy forma- 
tion and evolution. In particular, three landmarks have been 
pivotal in the development of theoretical studies of structure 
formation: a) the abundance of MW galaxy satellites that 
motivated one of the strongest points of tension with the 
now-standard A Cold Dark Matter (ACDM) paradigm of 



structure formation ( Klypin et al.|1999 Moore et al.|1999 l, 
b) the spatial distribution of the same satellites which trig- 
gered discussions on how unique the host dark matter halo 
of the MW is ( Metz et al.||2009 1 and c) the measurements 
of the tidal debris of disrupted merging galaxies around the 
MW and M31 galaxy, confirming the hierarchical nature of 
galaxy evolution, one of the fundamental characteristics of 
ACDM ( McConnachie et aL]|2009 K However, inferring gen- 
eral conclusions on galaxy evolution based on observations 
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of these two galaxies requires an assessment on how biased 
the properties of the MW and M31 are with respect to a 
given control population. 

In the framework of ACDM, the study of the MW and 
M31 starts by modeling their individual host dark matter 
halos, assuming that their simulated formation histories are 
"typical", or at least compatible with the assembly of the 
real Local Group ( de Rossi et al.|2009| Boylan-Kolchin et al. 
20101. The basic definition of a LG (in terms of the dark 



matter distribution) has two basic elements based on the 
state of the system today: (i) the estimated masses of the 
dark matter halos corresponding to the MW and M31 (see 
for instance Watkins et al. (20101 and references therein) 
and (ii) the isolation of these two halos from other massive 
structures ( Karachentsev et al.||2004 |. Two additional con- 
straints could be the separation and the relative velocity of 
the two halos ( Ribas et al.|20 05). However, the condition on 
the LG isolation admits a strict formulation, by requiring 
that the environment, in terms of the mass and position of 
the dominant galaxy clusters in the Local Universe, be as 
close as possible to the one inferred from observations. Such 
an additional condition imposes restrictions on the possible 
outcomes of structure formation on scales of the order of ~ 5 
Mpc. This is considered here as the meso-scale as opposed 
to the large (> 5 Mpc) or the small (< 1 Mpc) scales. 

The new feature in the analysis presented in this paper 
is the inclusion of such observational constraints around the 
LG environment in the initial conditions of the simulation. 
In a series of three simulations from such initial conditions, 
in a WMAP5 cosmology with a normalization ag = 0.817 
(Komatsu et al. 20091), we are able to define a sample of 



three LG dark matter halo pairs that form and evolve under 
specific conditions refiecting structure of the Local Universe. 
In addition we will take advantage of one of the largest cos- 
mological simulations carried out to date, the Bolshoi Sim- 



ulation (Klypin et al. 20101, to explore a larger sample of 



halos within the mass range of the LG, and calibrate possible 
cosmic variance effects. 

We analyze the constrained simulations with the pri- 
mary goal of quantifying the assembly histories of the LG 
halos. This is driven by two different motivations. One is to 
find out whether the simulated LGs, that are selected by 
dynamical considerations pertaining to their redshift zero 
structure, have mass aggregation histories (MAHs) that lead 
to the formation of disk galaxies like the MW and M31. The 
other is to find out whether such a MAH is dictated by by 
meso-scale environment of the LG, or whether a random se- 
lection of objects similar to the LG is likely to have a similar 
MAH. 

In Section [2J we describe our simulations and the 
method to re-construct the mass aggregation histories. In 
Section [3] we describe how we build the different control 
samples for our statistical analysis. In Section |4] we study 
the MAHs in the different samples and argue that the selec- 
tion by different isolation criteria does not induce a strong 
bias in the statistics describing the MAHs. In Section [S] we 
discuss the possible origin of these findings and comment on 
the connection with observations of the MW and M31. In 
Section |6] we summarize our conclusions. 



2 THE SIMULATIONS AND MASS 
AGGREGATION HISTORIES 

In this paper we make use of four cosmological N-body 
dark matter simulations. Three of them are part of the 
Constrained Local Universe Simulations (CLUES) project 
rl whose aim is to perform N-body cosmological simulations 
that reproduce the local large scale structure in the Uni- 
verse as accurately as current observations allow. The fourth 
simulation is the Bolshoi Simulation, which was performed 
from unconstrained initial conditions and spans a volume 
~ 60 times larger than each one of the CLUES simulations. 
In this section we will describe these simulations and the 
procedure we have used to construct the mass aggregation 
histories for the dark matter halos. 



2.1 The CLUES simulations 

First we describe the procedure employed to generate the 
constrained initial conditions. The observational constraints 
are the peculiar velocities drawn from the MA RK III|Willick 
et al. (19971, surface brightness fiuctuation Tonry et al 
( |2001 1 and the position and virial properties of nearby X- 
ray selected clusters of galaxies Reiprich & Bohringer ( 2002 1 



The Hoffman & Ribak ( 1991 1 algorithm is used to generate 



the initial conditions as constrained realizations of Gaus- 
sian random fields. These observational data sets impose 
constraints on the outcome of structure formation on scales 
larger than a few megaparsec. 

These constraints affect only the large and meso-scales 
of the initial conditions of the simulations, leaving the small 
scales essentially random. In particular, the presence of a lo- 
cal group with two dark matter halos roughly matching the 
masses, separation and relative velocities of the MW and 
M31 cannot be constrained. The strategy employed here is 
to construct an ensemble of 200 different realizations of the 
constrained initial conditions and simulate these with 256'^ 
particles on a box with side length 64/i~^Mpc using the Tree- 
PM MPI N-body code Gadget2 (ISpringelil2005l), and then 



scan these for appropriate LG-like objects within a search 
box centered on the actual position of the LG. Only three 
realizations are found to have such a LG object following the 
criteria detailed at the end of Sect.|3l It follows that the sim- 
ulations analyzed here obey two kinds o selection rules. By 
construction these are constrained simulations whose large 
and meso-scales are designed to mimic the local Universe. 
Then, post factum, the simulations that have the appropri- 
ate LGs are selected for further analysis. 

The selected simulations are then re-simulated at high 
resolution of 1024^ particles. The high resolution exten- 
sion of the low-resolution simulation is obtained by cre- 
ating an unconstrained realization at the desired resolu- 
tion, fast Fourier transforming it to fc-space and substitut- 
ing the unconstrained low k modes with the constrained 
ones. The resulting realization is made of unconstrained 
high k modes and constrained low k ones. The transitional 
scale happens around the length scale corresponding to the 
Nyquist frequency of the 256'^ mesh, Any = 2 x 64/256 
/i~^Mpc = 0.5/i^^Mpc . This corresponds to a mass scale 
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of M-My ^ 1.2 X 10^/i~^Mq, below which the structure for- 
mation can be considered as emerging primarily from the 
unconstrained k modes. 

The cosmological parameters in these high resolution 
simulations are consistent with a WMAP5 cosmology with 
a density Qm = 0.28, a cosmological constant Qa = 0.72, a 
dimensionless Hubble parameter h — 0.73, a spectral index 
of primordial density perturbations n = 0.96 and a nor- 
malization (78 — 0.817 ( _Komatsu et al. 20091. With these 
characteristics each particle has a mass rup — 1.89 x 10^ 
h-^UQ. 



2.2 The Bolshoi Simulation 



We have used as well the Bolshoi simulation (Klypin et al 



2010) to verify that the constrained simulation did not bias 
the halo samples and their MAHa~| The simulation was done 
in a cubic volume of 250 h~^Mpc on a side using 2048"' 
particles, leading to a particle mass of mp = 1.35 x 10* 
/i~^M0, roughly 10 times lower than the resolution in the 
CLUES simulations. 

We take from the Bolshoi simulation eight non- 
overlapping sub- volumes. Each sub- volume has a cubic size 
of 100 h~^Mpc on a side, corresponding to a comoving vol- 
ume comparable to the three CLUES simulations combined. 
The halo samples in the sub-volumes will be used to cali- 
brate the impact of cosmic variance on the different statistics 
we use to characterize the halo populations. 



2.3 Halo identification and merger tree 
construction 

In order to identify halos we use a FOF algorithm. We do 
not include any information of the substructure in each halo. 
All the analysis related to the mass aggregation history is 
done in terms of the host halos. In particular the mergers do 
not correspond to the fusion of an accreted sub-halo with a 
central dominant host halo, but instead correspond to the 
moment of two halos overlapping for the first time. 

The FOF algorithm has a linking length of b = 0.17 
times the mean inter particle separation. The mean over- 
density of objects found with this linking length at redshift 
z = Ois 680 ( [More et al.|2011[ ). We identify the halos for 80 
snapshots more or less equally spaced over the 13 Gyrs be- 
tween redshifts < z < 7. All the objects with 20 or more 
particles are kept in the halo catalogue and considered in 
the merger tree construction. This corresponds to a mini- 
mum halo mass of Mmin = 3.78 x 10^h~'^MQ. Within the 
CLUES simulations a Milky Way like dark matter halo of 
mass ~ 1.0 X lO^^/i~^M0 is resolved with ~ 5 x lO** parti- 
cles, in the Bolshoi simulation it is resolved with ~ 7 x lO"' 
particles. For the Bolshoi simulation we have used snapshots 
spaced by roughly 400Myr and followed the exact same pro- 
cedure to build the halo catalogues and the merger trees. 

Within the FOF analysis all FOF groups with 20 or 
more particles are identified. The merger tree construction 
is based on the comparison of the particles in FOF groups 
in two consecutive snapshots. Starting at z — for every 



FOF group in the catalog, Go, we find all the FOF groups 
in the previous snapshots that share at least thirteen parti- 
cles with Go and label them as tentative progenitors. Then, 
for each tentative progenitor, we find all the descendants 
sharing at least thirteen particles. Since the smallest FOF 
groups contain 20 particles at least 2/3 of the particles must 
be identified in tentative progenitors or descendants. Only 
the tentative progenitors that have as a main descendant 
the group Go are labeled as confirmed progenitors at that 
level. We iterate this procedure for each confirmed progen- 
itor, until the last available snapshot at high redshift. By 
construction, each halo in the tree can have only one de- 
scendant, but many progenitors. 

The mergers of FOF groups correspond to the time 
where the FOF radii of two halos overlap for the first time. 
The infall of the less massive halo into the host and the 
subsequent inspiral, disruption and fusion will be delayed 
with respect to the time of the FOF merger. Different the- 
oretical approximations and methodologies can predict the 
infall-fusion time-scale only as an order-of-magnitude esti- 



mate (Hopkins et al. 2010 1. The most used time-scale for this 



process is based on the Chandrasekhar dynamical friction 
formula, but improved estimates based on numerical sim- 



ulations ( Boylan-Kolchin et al 
yield 



2008 Hopkins et al. 20101 



tintall = 0.56 



(Mvir/M,, 



R^/ 



Vvir J ln(l -H Mvir/Meat) ' 



(1) 



where i?vir, K-ir and Mvir are the virial radius, velocity and 
mass of the host halo, Maat the mass of the future satellite at 
the moment of infall at 7?vir. A median initial circularity of 
the satellite orbit of 0.5 has been assumed. For mass ratios 

of Mvir/M,at = 10 



iinfall = 4.85 



-/Xvii 



(2) 



2.4 Local Group selection 

A LG in a constrained simulation consists of two main halos 
within a certain mass range, within a distance range and 
obeying some isolation conditions^ In addition it should re- 
side close to the relative position of the LG with respect to 
the Virgo cluster. Given the periodic boundary conditions 
of the simulations and the lack of treatment of the Zeldovich 
linear displacement in the reconstruction of the initial con- 
ditions, the large scale structure of the simulations is dis- 
placed by a few Megaparsecs among different realizations of 
the simulation. The most robust features of the constrained 
simulations are the Virgo cluster and the Local Supercluster. 
Their positions in the initial conditions are known, at 2 = 
their environment is searched for halos in the corresponding 
mass range to determine their present positions. These are 
used to fix the 'position' of the simulation in relation to the 
actual universe. In Table [T] we summarize the masses of the 
MW and M31 halos identified by the FOF halo finder in 
these three simulations. 



^ Halo catalogs for these simulations are available at 
Ihttp : //www . multidark . org/MultiDark/l 



•^ A quantitative description of these conditions is presented at 
the end of Section [sl 
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Figure 1. Halo distribution in the three CLUES and the Bolshoi (lower right) simulations. Only halos more massive than Mf^ > 
2 X lO^^/i^^MQhave been included. The radius of each circle corresponds to the radius defined by the FOF algorithm which is calculated 
to be the radius of a sphere with an equivalent volume as the FOF group. The dashed circle marks a 5h~^TsApc environment centered in 
the most massive halo of the LG. The solid thick circle shows the projected position of the halo identified with the Virgo cluster. The 
cut is 25/i~^Mpc thick and is centered at the LG position. 



Table 1. Properties of the MW-M31 pairs. Column 1: parent simulation. Column 2: Halo name (either MW or M31); column 3: FOF 
mass; column 4: last major merger time ; column 5: formation time; column 6: assembly time; column 7: matter overdensity calculated 
with in a sphere of 5/i~^Mpc . All times are look-back-times. 



Simulation 


Halo Name 


FOF Mass 


TM 


rp 


TA 


55 + 1 






[10 


^2ft-lM0] 


[Gyr] 


[Gyr] 


[Gyr] 




CLUES- 1 


M31 




1.39 


11.0 


11.0 


11.5 


0.72 


CLUES-1 


MW 




0.99 


10.0 


9.3 


9.7 


0.69 


CLUES-2 


M31 




0.98 


12.0 


10.0 


10.4 


0.78 


CLUES-2 


MW 




0.77 


11.3 


11.0 


11.0 


0.87 


CLUES-3 


M31 




1.45 


11.0 


10.6 


11.0 


0.75 


CLUES-3 


MW 




1.11 


9.8 


9.8 


11.0 


0.80 


Average 






1.15 


10.9 


10.3 


10.8 


0.76 


Standard Deviation 






0.23 


0.8 


0.6 


0.6 


0.05 
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FigurelTlshows the large scale structure of the three con- 
strained realizations centered on the position of the LG in 
each box in a slice 25/i~^Mpc thick. In the three CLUES sim- 
ulations shown in Fig. [Tithe projected position of the Virgo 
cluster is shown by a thick circle. The fourth panel in the 
same figure shows a cut of the same geometrical character- 
istics from the Bolshoi simulation, centered on one LG-like 
object. 



2.5 Merger trees description 

For each merger tree we define three different times to char- 
acterize the MAHs. Each time has direct connection with 
the expected properties of the baryonic component in the 
halo. The times, measured as look-back time in Gyr, are: 




• Last major merger time (tm): defined as the time 
when the last FOE halo interaction with ratio 1:10 starts. 
This limit is considered to be the mass ratio below which 
the merger contribution to the bulges can be estimated to 
be < 5%-10% ( [Hopkins et al.||2010| . Strictly speaking, as 
we do not follow sub-structure in the simulation, this event 
corresponds to the time when the merger fell into the larger 
halo and for the first time became a sub-halo. One can use 
Eq. ([2| to estimate the infall time-scale of the satellite to 
the center of the host. 

• Formation time (tf): marks the time when the main 
branch in the tree reached half of the halo mass at 2: = 0. 
This marks the epoch when approximately half of the total 
baryonic content in the halo could be already in place in a 
virialized object. 

• Assembly time {ta): defined as the time when the 
mass in progenitors more massive than M/ = 10^°/i~^MQis 
half of the halo mass at 2; = 0. This time is related to the 
epoch of stellar component assembly, as the total stellar 
mass depends on the integrated history of all progenitors 
(Neisteiii et al.||2006| [Li et al.||2008|). The exact value of 



TA is dependent on M/, the specific value selected in this 
work was chosen to allow the comparison of assembly times 
against the results of the Bolshoi simulation which has a 
lower mass resolution than the CLUES volumes. 

In Table [l] we summarize the values of these three dif- 
ferent times for the three pairs of MW-M31 halos. In Fig- 
ure [2| we show the median mass aggregation history in the 
main branch as a function of redshift for halos in the mass 
range5.OxlO"/i"^M0 < Mh < 5.0x lO^^/i"^M0. Following 



Boylan-Kolchin et al. (20101 we fit the MAH by a function 



of the kind 



M{z) = Mo(l + zf exp(-a(^rT^ - 1)), 



(3) 



with Q = 4.5 and /3 = 2.24. These values provide a good fit 
within 2.3% for 2 < 7. In the same Figure we overplot the 
main branch growth for the six halos in the three simulated 
LGs. The MAHs of these halos are systematically located 
above the mean, an indicator of early matter assembly with 
respect to the halos within the same mass range. 



Figure 2. Mass assembly histories of LG halos in the CLUES 
simulation as a function of redshift. The solid black line shows 
the median MAH for all halos in the CLUES simulations within 
the mass range 5.0 X IOI^/i-^Mq < M^ < 5.0 X IQI^/i-IMq, the 
dashed lines show the first and third quartiles. Also plotted as 
colour lines are the MAHs for the MW (dotted) and M31 halos 
(continuous) in the three constrained simulations. The assembly 
history for the LG halos is systematically located over the median 
values as sign of early assembly with respect to all halos in the 
same mass range. 



3 SELECTION OF LOCAL GROUPS AND 
CONTROL SAMPLES 

Four different samples of halos are constructed here, in a 
nested hierarchy in which the first sample contains the sec- 
ond which contains the third. The fourth sample is the one 
that includes the 3 LGs. These are to be used to study how 
the various criteria employed in constructing the samples 
affect the MAH of its members. The first three samples are 
constructed also from the Bolshoi simulation, and are used 
to look for possible biases in the constrained simulations. 

The first sample we define consists of all hal os in the 
mass range 5 x lO"/i"^M0< Mk <5 x 10^^ /t~^M0(l w"atkins 
&tal. 2010). We refer to this set as the Individuals halo 



sample. 

The second is a sample of halo pairs. Two halos, Ha and 
Hb , from the Individuals sample are considered a pair if and 
only if halo Hb is the closest halo to Ha and vice- versa. Fur- 
thermore, with respect to each halo in the pair there cannot 
be any halo more massive than 5.0 x 10^^/i~^M©closer than 
its companion ( [Karachentsev et al.||2004| ). We do not apply 
any further dynamical restrictions. For instance an element 
in this sample may be a pair of halos that are infalling into 
a cluster and are coincidentally close to each other. We refer 
to this set as the Pairs sample. 

The third is a sample of isolated pairs. We construct it 
by imposing additional conditions on each member of the 
previous sample. These conditions are defined to obtain a 
LG-like halo pair according to a series of requirements that 
follow the lines of [Governato et al.| ( [1997[ ), [Maccio et al.| 
( 2005[) [Martinez- Vaquero et al.[([2007[). We will refer to this 
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Table 2. Names and description of the four samples used to quantify the formation history of the LG halos. The three first samples 
are constructed both from the CLUES and Bolshoi simulations. By definition the LG sample can only be constructed from the CLUES 
simulations. The size refers to the total number of objects in the corresponding volume (individual halos or pairs). 



Name 



Description 



Size (CLUES) Size (Bolshoi) 



Individuals All the distinct halos in the mass range 5.0 X 10^^h~^MQ- 5.0 X lQ^'^h~^M(i 

Pairs All the pairs of halos constructed from the Individuals sample. 

Isolated Pairs Subset from the Pairs sample following some isolation criteria (see iJ3| 

LG The three pairs of LG halos from the constrained simulations. 
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Figure 3. Fraction of halos with merger histories described by a MAHs with rjv/, Tp and T4 larger than a given value. The lines 
represent different samples. The sample of Individuals (dashed) and Pairs (thick continuous lines) from the CLUES simulations and the 
Pairs extracted from eight sub- volumes in the Bolshoi simulation (thin continuous lines). 



sample as the Isolated Pairs sample. The conditions are the 
following: z 

a) The distance b etween the ce nter of the halos is smaller 
than 0.7 ft"^Mpc ( jRibas et al.|2005J . 

b) The relative radial velocity of the two halos is negative. 

c) There must not be objects more massive than either of 
the LG halos within a radius of 2/i~^Mpc from each object 
I TikhonovirKlypin|2()09| . 

d) There must not be a halo of mass > 5.0 x 
10^^/i~^Mq within a radius of 5/i~^Mpc with respect to each 
halo center ( Karachentsev et al. [20041. 



The final fourth sample contains the three objects that 
fulfil the criteria of third sample, and are located at about 
10/i~^Mpc 'south' of the Virgo cluster in the Supergalactic 
Plane. This sample is referred to as LG. 

We build these three samples both from the CLUES 
and Bolshoi simulations. A short summary description of 
each sample is contained in Table [2] 



4 RESULTS 

The backbone of our analysis is the study of the MAH of 
halos in the mass range (5.0xl0" - 5.0x lO^^)/i"^M0. Our 
results must be described in the 6 dimensional parameter 
space, spanned by the three characteristic times of the two 
halos, dubbed as MW and M31. The distribution of tm, 
tf and ta of the three different samples is studied in j ]4.1| 
and the possible dependence of these distributions on the 
ambient density around the LGs and the mass ratio of the 
MW and M31 members of the LGs in ■^n\ 



4.1 Mass accretion history of the different 
samples 

Figure [3] presents distribution of ta/, tf and ta for the Indi- 
viduals and Pairs samples of both the CLUES and Bolshoi 
simulations. The distribution with respect to the MW and 
M31 are virtually indistinguishable, and the curves present 
both halos. We calibrate the effect of cosmic variance with 
the (100ft~^Mpc)^ volumes extracted from the Bolshoi sim- 
ulation. The results are overplotted as thin magenta lines. 
The distribution of tm and tf are well within the scatter of 
the sub- volumes, while the ta is somewhat out of the range. 

We conclude that with respect to the MAH, the con- 
strained simulations are essentially unbiased with respect 
to the unconstrained one. The interesting fact that emerges 
here is the halos in the Individuals and Pairs sample share 
the same MAH, as expressed by the three times described 
here. 

Figure[4]presents the main results of the paper. It shows 
the distribution of the three times for the different sample 
of pairs of halos. The left column is made of three grey scale 
maps describing the number of objects in the Pairs sample 
in the subspace of (rf ^Srf ^), where X = M,F,A. The 
shades represent the number of pairs around a given region 
of parameter space calculated from the Pairs samples in the 
Bolshoi simulation. The three different LG pairs are over- 
plotted as stars. 

The right hand side column of Figure |4] shows the inte- 
grated relative distribution of the halos in the three different 
times of the Pairs, Isolated Pairs and LGs samples. For the 
LGs this is further separated for the MW and M31 halos. 
The distribution of the Isolated Pairs of the Bolshoi sub- 
volume is presented as well. 
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Figure 4. Left column. Joint distributions of tlirce different times (last major merger, formation and assembly) describing the mass 
aggregation histories. Each point in the plane represents a pair MW-M31 with histories described by the time values at that point. 
Levels in shading coding indicate the number of halo pairs in the Bolshoi simulations in that parameter range. Dark regions represent 
a high number of pairs. An absolute scaling for this shading can be obtained from the plots presented in the Right Column. The stars 
mark the location of the three LG pairs, each one coming from one of the constrained simulations. Right column Integrated probability 
of these three different MAH times. The continuous black lines represent the results for the Pairs sample in the CLUES simulations. 
The Isolated Pairs sample from CLUES is represented by the thick dashed lines. The results from the Isolated Pairs samples in eight 
sub-volumes of the Bolshoi simulation are represented by the thin continuous grey lines. The thick continuous lines represent the results 
for the LG sample. The distributions from the Pairs and Isolated Pairs control samples are basically indistinguishable. In other words, 
detailed selection criteria for halo pairs, based on isolation only, do not narrow down significantly the range of dark matter halo assembly 
properties. 
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Two important comments can be made based on Fig- 
ure HI First, we see that the times in the LGs sample are 
confined to a narrow range compared to the broad Pairs 
sample. The merger, formation and assembly times in this 
sample are confined within the range 9.5 — 12 Gyr. Second, 
from the integrated distribution, we infer that the Pairs and 
Isolated Pairs samples are virtually indistinguishable. This 



implies that the commonly used isolation criteria (Gover- 



nato et al.|1997||Macci6 et al.|2005l [Martinez- Vaquero et al 



2007) do not automatically produce the narrow parameter 



space occupied by the LG pairs. 



4.2 The influence of the Local Matter Density 
and the Mass Ratio 

The Pairs and Isolated Pairs samples are selected based on 
isolation and dynamics. The similarity of the distribution of 
the different MAH times of the the different samples moti- 
vates us to look for the possible dependence of these distri- 
butions on some other characteristics of the three LGs. In 
particular, the three LGs are found to share the two follow- 
ing properties: the mass ratio between the two halos and the 
matter over-density in a sphere of 5/i~^Mpc radiusrl noted 
as &C,. The values for the halo masses in the pairs and the 
local over-densities are listed in Table [l] together with the 
assembly, formation and last major merger times. A series 
of sub-samples of the halos sample are constructed by re- 
quiring that the masses and mass ratios between the pairs 
are bounded by the LG limits or the values 5,5. These sub- 
sampling do not bias the LG-like objects towards the region 
of parameter space defined by the LG sample. 



5 DISCUSSION 

Three basic facts emerge from the results presented in the 
previous section : a) the three LGs share a common for- 
mation history, b) this formation history is quiet out to at 
~ (10— 12)Gyr and c) none of the selection rules applied here 
to the pairs of halos have defined a sample of objects with 
MAH similar to that of the three LGs. In what follows, we 
discuss the possible origin and the predictable consequences 
of these facts. 



5.1 On the Common Formation History 

Naively, one might hypothesise that the fact that all three 
CLUES LGs have a common MAH, as defined here, is consis- 
tent with being drawn at random from the sample of pairs, 
i.e. the range of properties spanned by three random halo 
pairs can be naturally expected to be narrow. This is the 
null hypothesis we test now. 

What is the probability that 3 randomly selected pairs 
have MAHs within the range of properties found for the LG? 
We compute this probability based on the fraction of halos 
in the pair samples that share the LG formation properties. 

We define first the minimal subspace that contains the 



* ^5 has been calculated from the total mass in halos more mas- 
sive than 1 X IO^^/i^^Mq contained within a sphere of radius 
5/i~^Mpc centered at the position of each halo. 



3 simulated LGs by providing lower bounds on the differ- 
ent times describing the MAHs. Table |3] lists the minimal 
last major merger, formation and assembly look-back times, 
where two options are taken to estimate the minima. The 
first defines the "two sigma" bound, namely the average 
value minus twice the standard deviation of each time of 
the 6 halos of the 3 LGs, the second takes the minimum 
value for each time. 

The table provides the fraction of halos in the Individu- 
als sample satisfying each one of the conditions tx > t-^"""^ 
independently and all of them simultaneously, where X — 
M, F and A and the subscript bound denotes the minimal 
bound of such time. We find that the fraction of Individu- 
als in the quiet MAH subspace is fi = 0.17(0.12) both in 
CLUES and Bolshoi for the first (second) minima option. 
If we consider now the halos either in the Pairs or Isolated 
Pairs samples, only a fraction of fp = 0.03(0.01) pairs are 
composed of halos that are both within the LG parameter 
space. To a good approximation, the pair fraction can be 
calculated as the individual fraction squared, fp ~ fi x fi. 
This is the expected result under the assumption that the 
assembly of the MW and M31 are independent. 

The probability of randomly selecting three random 
halo pairs and having them within the range of parame- 
ters defined by the LG can be calculated as plg ~ fp ~ 
2.7 X 10"^(1.0 X lO"'^). This small probability is a conse- 
quence of having found 3 halo pairs within a set of proper- 
ties shared by 0.17(0.12) of the total population of halos. If 
we consider pairs with a range of desired properties within 
shared by, say, 0.68 of the halos in the total population (the 
fraction within one standard deviation around the mean), 
the probability of finding three pairs inside that range would 
bepi_^(0.68 X 0.68)3 ^ q ;^ 

Comparing the results of the probabilities plg and 
Pi~cr, the null hypothesis can be safely rejected. It is highly 
unlikely that the three randomly selected pairs show a nar- 
row range of properties as in the case of the LG sample. 

Both the ab-initio and post-factum constraints imposed 
on the LG yield a, LG sample with very sinrilar MAHs. In 
the CLUES simulations only the large and mid-scales are 
effectively constrained by the data leaving the galactic and 
smaller scales effectively random. It follows that the MAH 
of objects similar to the LG is strongly affected by their 
environment. To what extent this is valid for DM halos in 
general remains an open question. 

5.2 On the Quietness of the Formation History 

We established in the previous sections that the MAHs are 
quiet out to ~ (10 — 12)Gyr, and that none of the selection 
rules applied here to the pairs of halos have defined a sample 
of objects with MAH similar to that of the three LGs. 

The last point is consistent with the results previous 
studies that have approached the same question of estimat- 
ing a possible bias of the LG with respect to a general halo 



population ( de Rossi et al. 2009 Boylan-Kolchin et al. 20101. 
These studies apply isolation criteria on scales of lh~ Mpc 
over halos in the mass range we study here, and find as well 
that no significant bias is introduced in the isolated halo 
population with respect to the parent halo population. 

The parameter subspace defined by the three LG cannot 
by explained either in terms of the isolation criteria listed 
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Table 3. Fraction of halos/pairs the different samples with times tm, Tp and ta located in the parameter space defined by the minima 
characteristic times of the LG halos in the constrained simulations. These minima from the LGs are defined for each tx in two different 
ways: 1) as the mean value minus two times the standard deviation (see TablelTJ and 2) as the minimum value of all realisations. These 
minima times are denoted t'^ and t'^ respectively and are presented in the first rows. In the following rows, the first column describes 
the name and origin of the sample. The three following columns indicate the fraction of the total population with a tx larger than the 
calculated r^ or t'^ (in parenthesis). In the case of pairs samples, we require the times for both halos to be above the threshold. The 
last column refers to the three different tx being simultaneously larger than the corresponding r^ i'^x)- 



"Two sigma" bound 


t'm [Gyr] 
9.3 


r'p [Gyr] 
9.0 


t'a [Gyr] 
9.6 




"Minima" bound 


r'l, [Gyr] 
9.8 


r'^ [Gyr] 
9.3 


r'i [Gyr] 
9.7 




Sample 


TM > r;, (r^',) 


rp > r'p {r'^) 


-A > A (r'i) 


TM,F,A > t'm,F,A (^M,F,a) 



CLUES Individuals 
CLUES Pairs 
CLUES Isolated Pairs 
Bolshoi Individuals 
Bolshoi Pairs 
Bolshoi Isolated Pairs 



0.24 (0.18) 
0.06 (0.03) 
0.06 (0.03) 
0.23 (0.19) 
0.05 (0.04) 
0.05 (0.03) 



0.29 (0.24) 
0.09 (0.06) 
0.08 (0.05) 
0.23 (0.23) 
0.10 (0.05) 
0.10 (0.06) 



0.85 (0.85) 
0.74 (0.74) 
0.70 (0.70) 
0.87 (0.87) 
0.76 (0.76) 
0.73 (0.73) 



0.17 (0.12) 
0.03 (0.01) 
0.05 (0.03) 
0.17 (0.12) 
0.03 (0.02) 
0.03 (0.01) 



at the end of !|3] or by adding constraints on the values of 
the local over-density on 5/i.^^Mpc scales. The properties of 
the dynamical environment, common to all the CLUES sim- 
ulations and provide the quiet formation history for a LG, 
remain to be found. Ideally, that result should be expressed 
in a suitable form to search for LG pairs in an unconstrained 
simulation. 

Is the observed Local Group biased in the same man- 
ner? We cannot provide the answer to that question with 
the simulations we present in this paper. Nonetheless, the 
theoretical predictions we show here for the dark matter as- 
sembly in the LG seem to be in agreement with the disk 
dominated morphology of MW and M31. 



5.3 The Connection with the Observed Local 
Group 

The most distinct feature of the MW and M31 is that both 
galaxies have a disk dominated morphology. It is often men- 
tioned that abundant mergers, which are presumed to de- 
stroy the disk and be source of morphological change, are 
expected on all mass scales in the hierarchical picture of 
galaxy formation of ACDM generating a possible contra- 
diction with the abundance of disk galaxies in the local Uni- 
verse and, in particular, with the fact that the MW and M31 



A distinct and well characterised feature of the MW is 
the thick disk. This disk component of the MW has been 



are disk galaxies ||Toth fc Ostriker|l992 Quinn et al.|l993 



Kazantzidis et al. 2008 \ 



Our results provide new theoretical evidence that the 
MW and M31 could be expected to be disk dominated galax- 
ies in ACDM. From the results presented here, we have found 
that the last merger started on average 11 Gyr ago. At these 
redshifts the mass of the MW host halo is l-4x lO^/i'^M©, 
its virial velocity is « 200 km/s and its virial radius ~ 0.1 
h~^Mpc . Using these quantities and Eq.Mwe estimate the 
final infall time for the satellite to be « 3.5 Gyr, reaching the 
center « 7.5 Gyr ago. This quiet history should favour the 
survival of a disk formed in the halo ( Guedes et al.|[2011 |. 
Although, detailed estimations on these matters might have 
to include the inflow of gas into the disk ( [Scannapieco et aT" 
20091). 



known for more than 25 years ( Gilmore &i Reid||1983 K The 
thick disk contains a population of stars with different kine- 
matics, spatial distribution, ages and chemical enrichment 
compared to the thin galactic disk. Although M31 seems to 
have a similar component (Collins et al. 20111, the obser- 



vational and theoretical work on the MW's thick disk has 
a long history, and its origin can therefore be discussed in 
greater detail. One of the possible formation scenarios for 
the MW thick disk is an in-situ formation during/after a gas 
rich merger ( Sales et al.||2009 |. The analysis of the orbital 
eccentricity of stars based on RAVE and SDSS data sup- 
ports the gas-rich merger mechanism (Dierickx et al.|[2010 



Wilson et al. 20111. In our results the last merger reaches 



the center « 7.5 Gyr ago, close to the look-back time of 
SGyr as required by the in-situ formation scenario. 



6 CONCLUSIONS 

We use constrained simulations of the local Universe to 
study the dark matter mass aggregation history (MAH) 
of the Local Group (LG). Two basic questions motivate 
this study: 1. To what extent the simulated LGs can ac- 
count for the observed structure of the MW and M31 galax- 
ies? Namely, if the disk dominated morphology implies that 
the MW and M31 halos had a quiet MAH over the last 
~ llGyrs, can simulations recover this recent quiet history? 
2. Does this quiet MAH arise from the intrinsic properties of 
the DM halos, or is it induced by environment within which 
the LG is embedded? Is the implied MAH of the LG trig- 
gered by the large and meso-scales, or is it induced by the 
small, i.e. galactic and sub-galactic, scales? 

The methodology adopted here is to use constrained 
simulations of the local Universe, designed to reproduce the 
large and meso-scales of the LG environment, and search for 
halos that resemble the actual LG. The identification of a 
pair of halos as a LG-like object is based on a set of isolation 
and dynamical criteria, all formulated by their redshift zero 
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structure, in complete ignorance of their formation history. 
A LG-Uke object that is found close to the actual position 
of the observed LG with respect to the large scale structure 
environment is defined here as a LG. By construction a con- 
strained simulation can have only one LG or none at all. 
Indeed, out of a suit of 200 constrained simulations only 3 
harbour a LG. Controlled samples of individual halos and 
pairs have been constructed as reference samples. The anal- 
ysis has been extended to the unconstrained Bolshoi simu- 
lation that is used here for an unbiased reference (Klypin 
et al.|2010| ). 



The construction of the identification of the 3 LGs is 
done independently of the MAH of the halos. Yet, the MW's 
and M31's halos of the 3 LGs all have a common quiet MAH, 
defined as having the last major merger, formation and as- 
sembly look-back time extending over ~ (10 — 12) Gyr. 
This quiet formation history of the simulated LGs can help 
to explain the disk dominated morphology of the MW and 
M31, adding evidence to the internal instability origin of the 
spheroidal component of the MW ( Shen et al.|2010 \ . Based 
on measurements of the eccentricity of orbits in the MW, it 



has been recently claimed (jSales et al.|2009 Dierickx et al 



2010 Wilson et al. 20111, that a rich merger taking place 



10.5 to 8 Gyr ago is a favoured mechanism explain the thick 
disk in the MWlBrook et al.l (120041). Our finding of a quiet 



MAH of the LG provides a suitable platform for such a pro- 
cess to take place. 



The LG halos are assumed here to be selected from 
FOF halos in the mass range 5 x 10"/i"^Mo < Mh < 5 x 
10^^/i"^Mq at 2 = 0. Between 12% and 17% of these halos 
are found to have a quiet MAH, depending on the detailed 
definition of the quiet parameter space. From this point of 
view the MW and M31 halos are not rare. However, how 
likely is a pair of halos to have such a quiet history, shared 
by both halos? Making the naive null assumption that the 
MAH of a halo is an intrinsic property of a halo independent 
of its environment then the fraction of pairs should be the 
product of the fractions for a single halo. Indeed, the Pairs 
sample drawn out of the Bolshoi simulation confirms this 
assertion, finding that between 1% to 3% of the pairs have 
as quiet an MAH as the LG systems do. The probability of 
having selecting 3 pairs randomly and finding them with a 
quiet MAH is on the order of ~ 10""". 

Next, we look for what dynamical or environmental 
property determines rethe MAH of a LG-like object. We find 
here that the mere pairing of the MW-like halos does not 
affect the MAH fiducial times. Imposing the isolation and 
dynamical constraints that define the Isolated Pairs sam- 
ple does not affect it either. This leaves us with an open 
question as to what determines the MAH of halo pairs sim- 
ilar to the LG. The one hint that we have is that all of the 
three LGs reside in the same large and meso-scale environ- 
ment. We speculate that the cosmic web plays a major role 
in shaping the MAH of LG-like objects, although it is not 
yet clear what mechanism is responsible. A larger sample of 
constrained LGs is needed to confirm and further explore 
the reasons behind this result. 
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